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SUMMARY 


Lightweight aircraft structures exposed to a high intensity noise environ- 
ment can fatigue fail prematurely if adequate consideration is not given to the 
problem. Design methods and design criteria for sonic fatigue prevention have 
been developed based on analytical and experimental techniques. Most of the 
analytical work was based upon small deflection or linear structural theory 
which did not agree with the experimental results. A large deflection geomet- 
rical nonlinearity was incorporated into the analysis methods for determining 
the structural response to high intensity noise. The Karman-Herrmann large 
deflection equations with a single-mode Galerkin approximation, and the method 
of equivalent linearization were used to predict mean-square amplitude, mean- 
square stresses, and nonlinear frequency at various acoustic loadings for 
rectangular panels. Both simply supported and clamped support conditions with 
immovable or movable inplane edges are considered. Comparisons with experimen- 
tal results are presented. 


INTRODUCTION 


Vibrations caused by acoustic pressure can frequently disturb the operating 
conditions of various instruments and systems, and sonic fatigue failures which 
occurred in aircraft structural components cause large maintenance and inspec- 
tion burdens. The development of sonic fatigue data and design techniques were 
initiated to prevent sonic fatigue failures. Design methods and design criteria 
for many types of aircraft structures have been developed under Air Force 
sponsorship and by the industry in the past twenty years. Reference 1 has a 
complete list of the reports describing these efforts. This research led to 
sonic fatigue design criteria and design charts which are widely used during 
the design of an aircraft. Although current analytical sonic fatigue design 
methods are essentially based on small deflection or linear structural theory 
(see ref. 1, page 209), many documented tests (refs. 2-6) on various aircraft 
panels have indicated that high noise levels in excess of 110 decibels (dB) 
produce nonlinear behavior with large amplitudes of one to two times the 
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panel thickness in such structural panels. The neglect of such large deflec- 
tion geometrical nonlinearity in analysis and design formulations has been 
identified as one of the major causes for poor agreement between experimental 
data and analytical results. The evidence of those researchers was summarized 
in reference 7, where a conprehensive review of existing analytical methods on 
random excitations of nonlinear systems was also given. 

In this paper , the Karman-Herrmann large deflection equations for rectan- 
gular plates (ref. 8) are employed. Using a single-mode Galerkin's approxima- 
tion, the dynamic equations reduce to a nonlinear differential equation with 
time as the independent variable. The method of equivalent linearization 
(refs. 9 - 11) is then applied to reduce the nonlinear equation to an equivalent 
linear one. Mean-square displacements, mean-square stresses, and nonlinear 
frequencies at various acoustic loadings are obtained for rectangular panels of 
different aspect ratios and damping factors. Both simply supported and clamped 
boundary conditions with immovable or movable inplane edges are considered. 
Comparisons with experimental results are also presented. 


SYMBOLS 


a,b 

A,B 

^1 / ^2 
D 

err 

E 

f 

F 

h 

H(a>) 

L 

m 

N 

N 

P 

q 

r 

s (CO) 
t 

U, V 

w 

x,y,z 

B 

B* 

t 

A 

v 

P 

a,i 


Panel length and width 

Panel dimension parameters, 2n/a and 2Tr/b 

Constants 

Bending rigidity 

Error of linearization 

Young's modulus 

Equivalent linear frequency in Hz 

Stress function 

Panel thickness 

Frequency response function 

Spectrum level 

Mass coefficient 

Membrane stress resultant 

Constant 

Pressure loading 

Generalized or modal displacement 
Aspect ratio, a/b 

Spectral density function of excitation pressure p(t) 
Time 

Displacement of midplane 
Transverse deflection 
Coordinates 

Nonlinearity coefficient 

Nondimensional nonlinearity coefficient 
Ratio of damping to critical damping 
Nondimensional frequency parameter 
Poisson ' s ratio 
Panel mass density 
Normal and shear stresses 
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0) Radian frequency 

Equivalent linear or nonlinear radian frequency 
Subscripts : 
b Bending 

c Complementary solution or critical 

m Membrane 

max Maximum 

o Linear 

p Particular solution 


FORMULATION AND SOLUTION PROCEDURE 


Governing Equations 

Assuming that the effect of both the inplane and rotatory inertia forces 
can be neglected, the dynamic von Karman equations of a rectangular isotropic 
plate undergoing moderately large deflections are (refs. 8, 12) : 


L(w,F) = DV 4 w + phw fi _^_ - h (F # w 'xx + F ' 


— 2F , w, 

xy xy 


tt 

) - p(t) = 0 


xx w 'yy 


(1) 


V 4 F = E (w 2 , - w, w, ) 

xy 'xx yy 


(2) 


where w is the transverse deflection of the plate, h is the panel thickness, p 
is the mass density of the panel material, D = Eh 2 /12 (1-v 2 ) is the flexural 
rigidity, E is Young's modulus, V is Poisson's ratio, p(t) is the exciting 
pressure, and a comma preceding a subscript (s) indicates partial differentia- 
tion (s) . The stress function F is defined by 


a = f, 

x 'yy 

a = F, 
y xx 


( 3 ) 


T = -F, 
xy xy 


where O x , O y, and are membrane stresses. 

Simply Supported Panels . For a rectangular plate simply supported along 
all four edges as shown in Figure 1, Chu and Herrmann (ref. 8) , and Lin (ref. 

13) have considered that if the fundamental mode is predominant, the motion of 
the panel can be represented adequately as 

w = q(t) h cos (7Tx/a) cos (7Ty/b) (4) 

where q(t) is a function of time only. The maximum value of q(t) coincides with 
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the maximum deflection w max divided by panel thickness h. The expression w 
satisfies the boundary conditions for simple supports 

w = w, + V w, = 0 f on x = ± a/2 

f xx yy (5) 

w=w, + v w, = 0, on y = ± b/2 
yy xx 

Substituting the expression for w in Eq. (2) and solving for a particular 
solution Fp yields 

F = - ~r q 2 h 2 Er 2 (cos + -L cos ( 5 ) 

p 3 2 a r 4 D 

where r = a/b. The complementary solution to equation (2) is taken in the form 
- v 2 - x 2 

F = N V + N (7) 

c x 2 y 2 

where the constants N x and Ny contribute to the membrane stresses Q x and 0y and 
are to be determined from the inplane boundary, immovable or movable, 
conditions . 

For the immovable edges case, the conditions of zero inplane normal dis- 
placement at all four edges are satisfied in an averaged manner as 


7^- dxdy = [- (F, - VF,^) - ^ w 2 , x ] dxdy, on x = ± a/2 

J J ■* ■* 

dydx = j (F, xx - VF » yy > _ ^ w 2 , y ] dydx, on y = ± b/2 (8) 


where u and v are inplane displacements. For the movable edges case, the edges 
are free to move as a rigid body with the average inplane stress equal to zero. 
The inplane boundary conditions are 


N = h F, . 

X J-b/2 yY 


u = constant 
f a/2 


F , a y = 0l , 


on x = ± a/2 


N = h 

y 


F, dx =0 1, 


on y = ± b/2 


v = constant 

where N x and Ny are membrane stress resultants per unit length in plate. By 
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making use of these inplane edge boundary conditions, equations (8) and (9) , it 
easily can be shown that for the immovable edges 


_ _ q2 h 2 E „2 (1 + w2) 

X 8a 2 Cl-V 2 ) 

— q 2 h 2 ETT 2 o 

N = — (r 2 + V) 


y 8a 2 fl-V 2 ) 
and for the movable edges 


( 10 ) 


N = N = 0 
x y 


( 11 ) 


the complete stress function is then given by F = + F c - 

With the assumed w given by equation (4) and stress function given by 
equations (6) and (7), equation (1) is satisfied by applying Galerkin's method 


a/2 
^ -a/2 


b/2 

-b/2 


L(w,F) w dxdy = 0 


(12) 


from which yields the modal equation of the form 


q + U? q + 3q 3 = - P — 
^ 0 ^ ^ m 


(13) 


and 


to 2 - \ 2 — E — , A 2 — 77 4 (1 + ±-) 2 
° ° phb 4 ° * 2 


m = 7T 2 ph 2 /16 


(14) 


e - e + - (e* + e*) 

pc pc 


phb ' 


with 


* 3 tt 4 2 

6 « (1 + r ) (1 - v ) 

P a 4 
4r 


= —■ (1 + 2vr 2 + r 4 ) 
c 2r 4 


where 0) Q is linear radian frequency, m is mass coefficient, and 3 is nonlinear- 
ity coefficient. The linear frequency X 0 , nonlinearity coefficients 3p and 
6*, and aspect ratio r are all nondimens ional parameters. 
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Clamped Panels, Yamaki (ref. 14) considered the predominant mode 


w = a<t) n (1 + 


27 rx. 2 ttv. 

cos ) (1 + cos * ) 

t> 


( 15 ) 


which satisfies the clamped support conditions 


w = w, x = 0, 


w = w, = 0, 

y 


on x = ± a/2 
on y = ± b/2 


(16) 


By introducing equation (15) in equation (2) and solving it, the particular 
stress function is 

1 2 2 2 1 1 
— q h Er [cos Ax + — cos By + — cos 2 Ax 


32 

2 


+ — cos Ax cos By + yr cos 2 By 

(1 + 16r 4 

+ ' (4 + r2)"2 cos2Ax cos B y + - ( Y vVr^y 2 , COS ** COS 2 By ^ (17) 


where A = 27T/a and B = 27T/b. The complement airy stress function is assumed as 
the form appearing in equation (7) . Upon enforcing the inplane edge conditions, 
equations (.8) and (9) , it can be shown that for the immovable edges 


N 


x 


3q 2 h^E7T^ r - 2 \ 

— fl + vr) 

32 a 2 (1-V 2 ) 


N 


y 


3.q, 2 h 2 E7r, 2 (r 2 + v) 

32 a 2 (1-V 2 ) 


(18) 


and for the movable edges 


N = N = 0 
x y 


(19) 


the complete stress function is given by F = F p + F c . Introducing these ex- 
pressions for w and F in equation (1) and applying Galerkin's procedure yields 
the equation 


•• 

q 


+ 


2 

03 

0 


q + Bq 3 


pit) 

m 


(13) 
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where 


and 


03 = A 


2 D 


0 ^ 4 ’ 
phb 


x 2 = ml (3 + 2r 2 +' 3r 4 ) 
° 9r 4 


m = 9 ph 2 /16 


B - B p + B c - US* ♦ Bp ~~~Q 

phb 


* 4 4 2 r- 1 1 
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2 2 4 
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2 2 
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2 2 “ 

2(1 + 4r ) 


* 3 tt 2 4 

3 = (1 + 2Vr + r ) 

C n 4 

2r 


( 20 ) 


( 21 ) 


Equation (13) represents the undamped, large-amplitude vibration of a rectan- 
gular panel with simply supported or damped edges. 

The methods commonly used for determining the damping coefficient are the 
bandwidth method in which half-power widths are measured at modal resonances, 
and the decay rate method in which the logarithmic decrement of decaying modal 
response traces is measured. The values of damping ratio £ range from 
0.005 to 0.05 for the common type of panel construction used in aircraft 
structures. Once the damping coefficient is determined from experiments or 
from existing data of similar construction, the modal equation, equation (13) , 
now reads 


2 n 3 p(t) 

q + 2 ^ 0 ) q + 0 ) q+3q = ( 22) 

^ o o m 

The method of equivalent linearization is then employed to determine an 
approximate root-mean-square (RMS) displacement from equation (22) . 


Method of Equivalent Linearization 

The basic idea of the equivalent linearization (refs. 9-11) is to re- 
place the original nonlinear equation, equation (22) , with an equation of the 
form 


•• 

q + 


2C(0 Q q 


+ £2 q + err 


(q) 


P(t) 

m 


(23) 


147 



where ft is an equivalent linear or nonlinear frequency, and err is the error 
of linearization. An equivalent linear equation is obtained by omitting this 
error term, then equation (23) is linear and it can be readily solved. The 
error of linearization is 

err = (03^ - fi 2 ) q + (3q 3 (24) 


which is the difference between equation (22) and equation (23) . The smaller 
that the error is, the smaller the error in neglecting it, and the better 
approximate solution to equation (22) will be obtained. To this end, the 
equivalent linear frequency square ft 2 in the linearized equation is chosen in 
such a way that the mean-square error err 2 is minimized, that is 


9 (err 2 ) 

JW) 


0 


(25) 


If the acoustic pressure excitation p(t) is stationary Gaussian and ergodic, 
then the response q computed from the linearized equation, equation (23) , must 
also be Gaussian. Substituting equation (24) into equation (25) yields (refs. 
9, 13) 


ft 2 = (/ + 3$q 2 (26) 

where q 2 is the maximum mean-square deflection of the panel. Dividing both 
sides of equation (26) by D/phb^ yields 

2 2 * 2 

A = A q + 3g q (27) 

2 

where A is a nondimensional equivalent linear or nonlinear frequency 
parameter. 

An approximate solution of equation (23) is obtained by dropping the error 
term; the mean-square response of amplitude is 


q = 


S (03) H(CQ) I dOJ 


(28) 


where S (co) is the spectral density function of the excitation pressure p(t), 
and the frequency response function H(0)) is given by 

H (CO) = 1 (29) 

Itl(ft - 03 + 2i^03 Q 03) 

For lightly damped (£, < 0.05) structures, the response curves will be highly 
peaked at ft. The integration of equation (28) can be greatly simplified if the 
forcing spectral density function S (03) can be considered to be constant in the 
frequency band surrounding the nonlinear resonance peak ft, so that 
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C30) 


2 - tts m 

^ 2 2 
4m £0) R 

In practice, the spectral density function is generally given in terms of the 
frequency f in Hertz. To convert the previous result one must substitute 

ft = 27Tf (31) 

and S (ft) =: S(f)/27T 

into equation (30.) ; the mean^square peak deflection is simply 


32 S, 


2 = 


4 y . -v 2 

tt a 0 x 

32 S„ 


81 Z,X q V 


for simply supported panels 


for clamped panels 


(32). 


2 • 2 

The pressure spectral density function S(f)/2 7T has the units (Pa) /Hz or (psi) 
/Hz, and Sf is a nondimens ional forcing excitation spectral density parameter 
defined as 


S 


f 


sun 

2, 4, . ,.4,3/2 

p h (D/phb ) 


(33) 


The linear frequency parameters A Q in equations (32) are given in equation (14) 
and equation (.20) for simply supported and clamped panels, respectively, and 
the equivalent frequency parameters A 2 can be determined through equation (27) . 


Solution Procedure 


2 

The mean-square response q in equation (30) (or equation 32) is determined 
at the equivalent linear frequency ft (or A) which is in turn related to q^ 
through equation ('26) (or equation 27) . To determine the mean-square deflec- 
tion, an iterativ e p rocedure is introduced. One can estimate the initial mean- 
square deflection q2 using linear frequency G0 o through equation (30) as 


2 


q 


o 


TTS (a ) 0 ) 

~2~ 3 

4m £0) o 


C34) 


This initial estimate of q£ is simply the mean-square response based on linear 
theory. This initial estimate of "q? can now be used to obtain refined estimate 
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2 2 2 2 

of £1 through equation (26), £1=0) + 3g q , then q, is obtained through 

1 • « \ 1 o o 1 

equation (30) as 


— 7TS (fix) 

< 3i — 2 o 

4hi cco^ 2 


( 35 ) 


As the iterative process converges on the n-th cycle, the relation 


2 


TFS (flp) 

4m 2 C0) o fi 2 n 



1 


(36) 


becomes satisfied*, In the numerical results presented in the following section, 
convergence is considered achieved whenever the difference of the RMS displace- 
ments satisfied the relation 



(37) 


Stress Response 

Once the RMS displacement is determined, the bending stresses on the sur- 
face of the panel can be determined from 

6D 


O (w, + Vw, ) 

xb ^2 xx yy 


6D , , x 

7 ^ (w, + Vw, ) 

yb 2 yy xx 


(38) 


From equations (3) and (38), and using equations (4), (6), (7) and (10), the 

expressions for the nondimensional stresses on the surface of a simply supported 
panel with immovable edges are given by 

2 

Ob 2 2 

X 2 = (a + a ) — — = [ — — + V) cos — cos q 

Eh 2 xb xm Eh 2 2 (1-v 2 ) r 2 a b 


, , TT 2iry, 2 . it (1 + Vr ) , 2 

+ (—3- cos -r*0q + [ 5 — ] q 

8r 8r (1-v 2 ) 


750 
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■ - ~ = (CT + O ) 
Eh 2 yb 


Eh 
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IT 


2(1-V 2 ) 


v_ 

2 J 


ttx Try t 

cos rr“' 

a b 


(1 + -~r-) cos cos *^] q 


2 2 2 
. ,7T 2ttx v 2 . r 7T (r +V) , 2 

+ (— cos — -) q + [ 5 —l q 

8r (1-V ) 


( 39 ) 


For movable inplane edges, the last term in equation (39) vanishes. Simi- 
larly, from equations (3) and (38), and using equations (7), (15), (17), and 

(18) , the expressions for the nondimens ion al tensile stresses on the surface 
of a clamped panel with immovable edges are 


a b 
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Eh 2 (1-V ) r 


TT 1 

— - — [— — cos Ax(l+cos By) + V(l+cos Ax) cos By] q 


2 2 , 

. tt r r 1 

+ — — [— cos By + 
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2 2 
(1+r ) 


cos Ax cos By + — cos 2 By 
4r 


2 x 2 
(4+r ) 


cos 2 Ax cos By + 


(l+4r 2 ) 


cos Ax cos 2By]q^ 


+ [ 


3TT 2 (l+vr 2 ) 


32r (1-V ) 


■] q 


(40) 


a b Tr 2 v 

y = — [ (l+cos Ax) cos By + — cos Ax(l+cos By)] q 

Eh 2 (1-V ) r 
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4 12 
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2 2 
(4+r ) 

3JT 2 (r 2 +V) 2 

+ [ 2 2 q 

32 r (1-V ) 


2 2 
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where A = 27T/a and B = 27T/b. For movable edges, the last term in equation (40) 
vanishes. 
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Examining equations (39) and (40), a general expression is obtained for 
the stress at any point in the structure as 

a = C^q + C 2 q 2 (41) 

where and C 2 are constants. The constants can be determined from material 
properties, dimensions of the panel, and the location and direction at which 
the stress is to be measured. The mean-square stress is then related to the 
mean-square modal amplitude in a general expression as 

~T 2 ~~2 2 ~2 2 

O = q + 3 C 2 (q ) (42) 

~ 2 " 

Once the mean-square deflection q is determined, equations (36) and (37), the 
mean-square stress can then be obtained from equation (42). 


RESULTS AND DISCUSSION 


Because of the complications in analysis of the many coupled modes, only 
one-mode approximation is used in the formulation. The assumption for funda- 
mental mode predominacy is admittedly overly simplified; the conditions under 
which this is a valid approximation remain to be investigated. However, a 
simple model sometimes helps to give basic understanding of the problem. 

Using the present formulation, response of nonlinear rectangular panels 
with all edges simply supported and all edges clamped subjected to broadband 
random acoustic excitation are studied. Both immovable and movable inplane 
edges are considered. In the results presented, the spectral density function 
of the excitation pressure S(f) is considered flat within a certain region 
near the equivalent linear frequency f and a value of Poisson's ratio of 0.3 
is used in all computations, unless otherwise mentioned. Mean-square ampli- 
tudes and mean-square nondimens ional stresses for panels of various aspect 
ratios and damping ratios are determined and presented in graphical form. 

These graphs can be used as guides for preliminary design of aircraft panels. 

The maximum mean-square deflection can be reasonably obtained from these 
figures; however, multiple modes had to be considered for accurate determina- 
tion of mean-square stresses. This has been demonstrated by Seide in 
reference 15 for a simple beam subjected to uniform pressure excitation and in 
reference 16 for large deflections of prestressed simply supported rectangular 
plates under static uniform pressure. Comparison with experiment is also given. 
It is demonstrated that the present formulation gives remarkable improvement in 
predicating RMS responses as compared with using the linear theory. 

Analytical Results 

Figure 2 shows the maximum mean-square nondimens ional deflection versus 
nondimensional spectral density parameter of acoustic pressure excitation for 
rectangular panels of aspect ratios r = 1, 2, and 4, and a damping ratio 0.02. 

It is clear from the figure that an increase of r will "close" the curve. 
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This occurs because as r increases, the panel becomes less stiff, and the 
mean-square deflection has to be finite. It can also be seen from the figure 
that the mean-square deflection of the movable inplane edges case is approxi- 
mately twice as that of the immovable edges. 

The maximum mean-square nondimensional stress (bending plus membrane 
stress, at the center of the panel and in the y-direction) is given in Figure 
3 as a function of excitation spectral density parameter for simply supported 
rectangular panels of various aspect ratios and a damping factor 0.02. 

Results showed that the difference of maximum mean-square stresses between 
immovable and movable edges is small as compared with the difference of mean- 
square deflections between the two edge conditions. 


Figure 4 shows the mean-square deflection versus forcing spectral density 
parameter for simply supported square panels of different damping ratios. The 
corresponding maximum mean-square stress (bending plus membrane stress, at the 
center of panel) is shown in Figure 5. As it can be seen from the figure that 
the precise determination of damping ratio from experiment is important, e.g., 
stress increases by 25-30 percent as t, is decreased from 0.015 to 0.01 (for 
between 5000 to 20000) . 


Plots of the equivalent linear or nonlinear frequency parameter A 
versus mean-square modal amplitude for simply supported rectangular panels of* 
aspect ratios r = 1, 2, and 4 are shown in Figure 6. The lowest value of A 2 
corresponds to the linear case. 


In Figure 7, the mean-square deflection is given as a function of excita- 
tion spectral density parameter for rectangular panels of aspect ratios 
r = 1, 2, and 4 and a damping ratio 0.02. The maximum mean-square deflection 
of the clamped panels is somewhat much less than that of the simply supported. 
The corresponding maximum mean-square nondimensional stress (bending plus 
membrane stress, in the y-direction and at the center of the long edge) versus 
spectral density parameter is shown in Figure 8. 

Figure 9 shows the mean-square modal amplitude versus spectral density 
parameter of excitation for a square panel of different damping ratios. In 
Figure 10, the equivalent linear frequency parameter is given as a function 
of mean-square deflection for clamped rectangular panels of aspect ratios 
r = 1, 2 , and 4. 


Comparison with Experimental Results 

The experimental measurements on skin-stringer panels exposed to random 
pressure loads reported in references 3 and 4 are used to demonstrate the 
improvement in predicting panel responses by using the present formulation. 
The structure was a skin-stringer, 3-bay panel as shown in Figure 11. The 
panels were constructed of 7075-T6 aluminum alloy. Details of the test 
facility, noise sources, test fixture, and test results are given in 
reference 3. The important properties of the panel are 
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Length 


a = 68.58 cm (27 in.) 


Width between the rivet lines b = 16.84 cm (6.63 in.) 

Thickness h = 0.81 mm (0.032 in.) 

Damping ratio £ = 0.0227 

Poisson's ratio V = 0.33 

3 6 

Young's modulus E = 66.19x10 MPa (9.6x10 psi) 

Weight density p = 7.164 kg/m 3 (0.1 lb/ in. 3 ) 

The tests were conducted with an overall sound pressure level (SPL) of 157 dB, 
with a range of +1-5 dB which corresponds to an average spectrum level of 
125.26 dB (see Table IV of ref. 3 or Table 8 of ref. 17). The central bay of 
the 3-bay test panels is simulated by a flat rectangular plate. The linear 
frequencies for both simply supported (equation (14)) and clamped (equation 
(20)) support conditions are calculated and shown in Table 1. Test measure- 
ments and finite element solution are also given for comparison. Table 1 also 
shows the equivalent linear or nonlinear frequencies at overall SPL 157 dB. 

Table 1. Frequency Comparison 

Natural Equivalent linear 

frequency f Q frequency 


Simply supported - Immovable edges 

71 

321 

- Movable edges 

71 

240 

Clamped - Immovable edges 

159 

311 

- Movable edges 

159 

264 

Finite element (ref. 4) 

155 

N/A 

Experiment (ref. 3) 

126 , 129 

N/A 


Frequency at high intensity noise level was not reported in reference 3. From 
the results shown in Table 1, it is clear that the central bay of the test 
panels did not respond to the acoustic excitation as though it were fully 
clamped on all four edges. This was also demonstrated in Figures 12 and 17 of 
reference 3 in the sense that the highest measured RMS strains did not occur 
at the center of the long edges. The central bay of the test panels actually 
behaved somewhat between fully simply supported and fully clamped support 
conditions . 

The acoustic pressure spectral density S(f) is related to the spectrum 
level L as 
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J 8.41 x 10 (L/1 ° " 18) (psi) 2 /Hz 

S(f) = < (47) 

I 4 x 10 (L/1 ° ~ 8) (dynes/cm 2 ) 2 /Hz 

A spatially uniform white noise pressure loading with spectral density of 
S(f) = 2.824 x 10~ 5 (psi)^/Hz (or nondimensional spectral density parameter 
S f = 5100), which corresponds to an average spectrum level L = 125.26 dB, is 
used in the computations. The RMS stresses (equation (42)) at the center of 
the long edges for simply supported (equation (39)) and clamped (equation (40)) 
boundary conditions are calculated and given in Table 2. 

Table 2. Stress Comparison 

(RMS stress, kpsi at overall SPL 157 dB) 



ft 

T 

/ 

IF 


3 

< 

Ay/ 

y 


Linear 

Nonlinear 

Linear 

Nonlinear 


Theory 

Theory 

Theory 

Theory 

Simply-Supported 

0.0 

0.58 (Im. ) 

0.0 

3 .28 (Im. ) 



0.17 (Movable) 


2.74 (Movable) 

Clamped 

2.17 

1.12 (Im. ) 

6.57 

3.84 (Im. ) 



1.32 (Movable) 


4.24 (Movable) 

Finite Element (ref. 4) 

2.4 

NA 

7.7 

NA 

Experiment (refs. 3, 4) 





Panel A 


0.63 


2.2 

Panel B 


0.94 


2.9 

Panel C 


0.78 


2.5 

Panel D 


1.1 


- 

Panel E 


0.84 


2.2 

Average A-E 


0.87 


2.5 


Table 3 shows the RMS deflections using the present formulation. The 
measured and finite element RMS stresses and RMS deflections in reference 4 
are also given in the tables for comparison. It demonstrates that a better 
correlation between theory and experiment can be achieved when large deflec- 
tion geometrical nonlinearity effect is included in the formulation. 
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Table 3. Deflection Comparison 


/ (W /h ) 2 

/J max 



Linear Theory 

Nonlinear Theory 

Simply Supported 

8.0 

1.8 

(Immovable) 



2.4 

(Movable) 

Clamped 

2.7 

1.4 

(Immovable) 



1.6 

(Movable) 

Finite Element (ref. 4) 

3.1 

NA 


Measured (refs. 3, 4) 

- 

2.0 



CONCLUDING REMARKS 


An analytical method for predicating response of rectangular nonlinear 
structural panels subjected to broadband random acoustic excitation is pre- 
sented. The formulation is based on the Karman- Herrmann large deflection 
plate equations, a single-mode Galerkin approximation, the equivalent linear- 
ization method, and an iterative procedure. Both simply supported and clamped 
support conditions with immovable or movable inplane edges are considered. 

Panel mean-square deflection, maximum mean-square stress, and equivalent linear 
frequency at given excitation pressure spectral density can be determined, and 
they are presented in graphical form. These graphs can be used as guides for 
preliminary design of aircraft panels under high noise environment. Results 
obtained agree well with the experiment. It is suggested that further research 
be carried out with special attention to employ multiple modes in the formula- 
tion for accurate determination of mean-square stresses, and additional test 
data on simple panels are needed for an adequate quantitative comparison. 
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Figure 3. Maximum mean-square stress versus spectral 
density parameter of excitation for simply 
supported panels , £ = 0-02. 
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Figure 4. Effect of damping on mean-square deflection for a 
simply supported square panel. 
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SECTION A-A 

NOTE: DIMENSIONS IN INCHES 

Figure 11. Skin-stringer panel (after Van der Heyde and 
Smith, ref. 3). 


163 



